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Cooper pairing instability in a Fermi liquid is well understood by the BCS theory, but pairing mechanism for 
doped Mott insulators still remains elusive. Previously it has been shown by density matrix renormalization 
group (DMRG) method that a single doped hole is always self-localized due to the quantum destructive 
interference of the phase string signs hidden in the f-/ladders. Here we report a DMRG investigation of hole 
binding in the same model, where a novel pairing-glue scheme beyond the BCS realm is discovered. 
Specifically, we show that, in addition to spin pairing due to superexchange interaction, the strong 
frustration of the phase string signs on the kinetic energy gets effectively removed by pairing the charges, 
which results in strong binding of two holes. By contrast, if the phase string signs are "switched off 
artificially, the pairing strength diminishes significantly even if the superexchange coupling remains the 
same. In the latter, unpaired holes behave like coherent quasiparticles with pairing drastically weakened, 
whose sole origin may be attributed to the resonating-valence-bond (RVB) pairing of spins. Such non-BCS 
pairing mechanism is therefore beyond the RVB picture and may shed important light on the high- T c 
cuprate superconductors. 

The Cooper pairing is the hallmark of superconductivity in both the conventional and the unconventional 
superconductors as evidenced by experiments like, e.g., the flux quantization. In the BCS theory, two 
electrons injected into a Fermi liquid can always form a Cooper pair under an arbitrarily weak attractive 
interaction. Ever since the discovery of the high-T c cuprates, a great effort has been devoted to finding the 
responsible pairing glue, which is widely attributed to the superexchange interaction 16 . 

However, at a deeper level, the BCS theory as a suitable description of the Cooper pairing has been seriously 
challenged in the cuprate 7 . For instance, in a Mott insulator, the strong on-site Coulomb repulsion will cause the 
charge being stripped off the electrons, while their spins form the singlet RVB pairing in the RVB theory 8 " 10 . 
Furthermore, in a doped Mott insulator, the original Fermion sign structure for a non-interacting Fermi gas is 
replaced by a much sparse sign structure (phase strings) 11 " 14 as precisely identified in the {-/model at arbitrary 
dimensions. 

Without the integrity of individual electrons, to understand the nature of Cooper pairing in the doped cuprates, 
one has to go beyond the BCS scheme of simply identifying the pairing glue. Here, the behavior of the unpaired 
single-particle excitation has to be examined simultaneously. Recently, the loss of quasiparticle coherence in a 
Mott insulator has been studied for the single hole doped t-J square ladders by a large-scale DMRG simulation 15 . 
Due to the destructive quantum interference effect of phase strings, a novel charge localization purely of strong 
correlation origin has been unveiled, which is independent of whether the underlying spin correlation is quasi- 
long-ranged (in odd-leg ladders) or short-ranged (in even-leg ladders). It is also found that a coherent Bloch 
quasiparticle behavior can be recovered once the phase string effect is artificially turned off in the kinetic term of 
the t-J model 15 . 

Experimentally, the absence of a coherent quasiparticle excitation has been clearly observed by the angle- 
resolved photoemission spectroscopy (ARPES) in lightly doped cuprates such as Ca 2 - x Na x Cu0 2 Cl 2 16,17 as well as 
the underdoped YBa 2 Cu 3 0 6 + x 18 . As a matter of fact, the single-particle excitation is generally frustrated in a 
normal state of the cuprate, from the antiferromagnetically (AF) ordered phase to the pseudogap regime, and to 
the strange-metal phase at the optimal doping 1 ". The transport experiment has also universally shown the 
localization of charge carriers in the underdoped regime before superconductivity sets in 20 . 

In this paper, we further investigate the hole pairing in the t-J ladders using the DMRG algorithm 21 . A 
significantly strong binding energy is indeed found for two holes injected into a short-range-ordered (even- 
leg) spin ladder, but the pairing strength becomes vanishingly small for two holes doped into a gapless (odd-leg) 
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Figure 1 | The binding energy Ef, for two holes in a f-/ladder is obtained by the DMRG method, with the leg number ranging from 1 in (a), 3 in (b), 
and 2 in (c), 4 in (d). For even-leg ladders with finite spin gaps, the binding strength is substantially larger than the odd-leg ladders without a spin gap. In 
the large length limit, the extrapolated binding energy almost vanishes for the odd-leg ladders, while remains big for the even-leg ladders until at 
much larger f //2> 1 ■ It is interesting to note that for the even-leg ladders, the binding strength is strong not only in the singlet channel, but also in the triplet 
(SL = 1) channel. 



spin ladder. The origin of pairing here supports an RVB picture 
rather than the conventional BCS picture of exchanging the magnetic 
fluctuation. Namely, the hole pairing naturally arises from doping a 
short-ranged spin liquid state. Further surprisingly, the short-range 
spin correlation is only a necessary condition, but not a sufficient 
one. By turning off the phase string effect in the hopping term with- 
out changing the spin correlation, the coherent Bloch state is restored 
for the unpaired holes. But simultaneously the strong pairing disap- 
pears as well. It unveils a novel Cooper pairing mechanism, which 
works as an intrinsic combination of the spin RVB pairing with the 
charge pairing that removes the kinetic energy frustration of phase 
strings. Moreover, we verify that bound pairs of holes are generally 
repulsive to each other without further forming droplet. Our model 
study points to a non-BCS route to achieve high- 1",, in a doped Mott 
insulator, i.e., making the AF correlation a short-ranged one and, at 
the same time, frustrating the kinetic energy as much as possible by 
phase strings in the normal state. 

Results 

Model Hamiltonians. As a large- U Hubbard model with the 
hopping integral fS>/, where the superexchange coupling / = 4P/ 
U, the t-J Hamiltonian is given by H t .j = H t + Hj with 

1 



where (if) stands for the nearest neighbors (NN). c ia is the electron 
creation operator at site i, S ; and n, are the spin and number 
operators, respectively. The Hilbert space is constrained by the no- 
double-occupancy condition, i.e., n t £ 1. Our study focuses on the 
ladders on square bipartite lattices of N = N x X N r where N x and N y 
are the site numbers in the x and y directions, respectively. For the 
present DMRG simulation, we set / = 1 as the unit of energy. 



Binding energy. Define the binding energy by 
Eb = E 2 +E 0 — 2Ei , 



(2) 
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where E 2 and Ei are the ground-state energies of the two-hole and 
one-hole states, respectively, while £ 0 denotes the ground-state 
energy at half filling. If two holes form a bound state, then Ef, is 
negative as E 2 — E 0 < 2{Ei — E 0 ). Otherwise, the binding energy 
Eb should vanish in the thermodynamic limit with E 2 — E 0 = 2{Ei — 
E 0 ) for two independent holes. 

The binding energy Ey is studied by the DMRG method in the two- 
hole doped t-J ladders, with the leg number ranging from N y = 1 to 4 
as shown in Fig. 1. The binding strength is substantially large for an 
even-leg ladders [i.e., N y = 2 and 4 in Figs. 1 (c) and (d), respectively] , 
whereas Ey, diminishes quickly for an odd-leg ladder [i.e., N y = 1 and 
3 in Figs. 1 (a) and (b), respectively] with the increase of the ladder 
length N x . Indeed, for the odd-leg ladder cases, E b can be extrapolated 
to a vanishingly small value in the thermodynamic limit according to 
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Figure 2 | The contributions of the hoping and superexchange terms to the binding energy are separated by E b = {H t ) b + (H^ b : (a) 2-leg ladder in the 
singlet (S z tot = 0) channel, (b) 2-leg ladder in the triplet (S z tBt = 1) channel. 



a finite size scaling using second-order polynomials of l/N (see 
Supplementary Materials). 

Note that the main distinction between the even- and odd-leg 
ladders is well known for the undoped case: there is a robust spin 
gap in the even-leg ladders with exponential-decay spin correlations 
at a length scale about 2 to 3 lattice spacing for N y = 2, 4, but the spin 
excitation is gapless for the odd-leg ladders with quasi-long-range 
spin correlations 22,23 . It suggests that a short-range AF correlation in 
the spin background should be a necessary condition for a meaning- 
ful pairing strength as found in Fig. 1. Such an even-odd distinction 
in pairing has already been seen in the previous DMRG work at 
smaller lattice sizes 3 . As a matter of fact, it supports an RVB picture 810 
of pairing: two doped holes gain a binding energy by removing an 
RVB pair from the spin background. As for a weaker but still sub- 
stantial binding energy for the triplet pairing, shown in Figs. 1 (c) and 
(d), it means that the lowest energy of a triplet spin excitation still lies 
below two free spinon excitations in the even-leg ladders 24,25 . By 
contrast, in the odd-leg ladders, the RVB pairing is long- ranged with 
gapless free spinon excitations, such that the binding energy is van- 
ishingly small in the thermodynamic limit. 

One can further distinguish the contributions from the kinetic and 
superexchange terms in the binding energy. For this purpose, we 
redefine E b = {H t ) b + {H r ) b , where {H t ) b = {H t ) 2 - 2(H ( > 1 and 
(H,) b = (Hj) 2 + <H 7 > 0 - 2<ff J ) 1 . Here (H t ) 2 «H 7 > 2 ) and {H t \ 
((Hj)^ represent the kinetic energy (superexchange energy) of the 
two-hole and one-hole doped systems, respectively, and (Hj) 0 the 
superexchange energy at half filling. In Figs. 2 (a) and 2 (b), the 



separated contributions for the two-leg ladder are shown in the 
singlet and triplet channels, respectively. Clearly, the superexchange 
interaction serves as the driving force for the hole binding with (Hj) b 
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Figure 3 | The binding energy of the two pairs in the four hole case, 
defined as £j|-h°ies _ (jg, _ 2(E 2 -Eq). It shows that two hole pairs are 
repulsive for all the finite systems, indicating that the binding exists only 
between two holes. 
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Figure 4 | (a) By inserting a flux <1> into a closed ribbon made of a two-leg ladder connected along the long chain direction (the inset), the energy 
difference between <1> = 0 and<l> # 0: A£^ holes (O) = £ 2 (<I>) —£2(0), is calculated at f = 3/. Here<l> = n corresponds to the flux quantization condition, and 
A£^ holes oscillates strongly and decays exponentially with a length scale X ~ 3, indicating the frustration of the phase string effect and the pairing size. The 
binding energy Et, as a function of N x is also presented, (b) AEq oles at a non-quantized $ = k/2 exhibits a power-law decay at large N x (the inset), 
indicating that the centre-of-mass motion of the hole pair behaves like a phase-string-free coherent object, which is not shown in (a) because its 
contribution at <S> = n is the same as O = 0 (see the text), (c) For the single-hole case, A£j," holes (<l>) =£i(0) —£^0) also exhibits an oscillation with an 
envelop of exponential decay, indicating 15 the self-localization of the hole with { = 14.5 due to the phase string effect. 



< 0 whereas (H t )i, > 0. The hole binding eventually abruptly 
vanishes at t/J > 10 (singlet pairing) and t/J > 5 (triplet pairing) 
for the two-leg ladder in Fig. 2. Such abruptness is actually consist- 
ent with an RVB picture rather than a conventional pairing mech- 
anism by exchanging magnetic fluctuations. In the latter, a smooth 
crossover to the disappearance of the Cooper pair is usually 



expected with the increase of t/J. Once the total binding energy 
equals to zero, one has (H f )j, = (if;);, = 0 such that the two doped 
holes behave independently. Note that at t < /, the binding energy 
£(, is found to quickly reduce as the holes tend to stay at the two 
sides of the open boundary in the a DMRG calculation. But the 
boundary effect and phase separation at t/J < 1 are no longer 
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Figure 5 | (a) The motion of a single hole (open circle) generally leaves a 
sequence of signs (i.e., the phase string) on its path, besides an unpaired 
spin (red arrow) in the spin background, (b) The binding of two holes can 
remove the unpaired spins to gain superexchange energy, but at the same 
time the frustration due to phase strings is compensated as well to gain the 
kinetic energy. Here a phase string is a product of the signs depending on 
the background spins exchanged with the hole during its hopping, and the 
short-range pairings of them effectively eliminate the destructive 
interference effect of phase strings. 

important once the kinetic energy becomes dominant over the 
superexchange energy at t > /, as illustrated in Figs. 1 and 2, which 
is our main focus in this work. 

Finally, we check that there is no formation of a "droplet" when 
more holes are added. Define the binding energy for two pairs of 
holes by ££- holes = (£ 4 -£o)-2(£ 2 -£o), where £ 4 is the ground- 
state energy of the four-hole state. As shown in Fig. 3 for the two- 
leg t-J ladder, the two hole pairs are actually repulsive to each other at 
any finite size, and do not form a 4-hole droplet in the thermodyn- 
amic limit. Since each Cooper pair is well formed in a spin gapped 
state (its coherence will be further examined below), a superconduct- 
ing condensation is naturally expected for a finite density of holes. 

Novel pairing mechanism. In the following, we show that the 
presence of a spin gap/short-range spin correlation (i.e., the RVB 
mechanism) is only a necessary but not a sufficient condition for 
the appearance of strong binding between the holes. Instead, a new 
pairing mechanism hidden in the kinetic energy term of the t-J model 
will play an essential role, which is of non-BCS type. 

To examine the nature of pairing, we connect the ladder along the 
x direction to make a close loop of circumference N x and then thread 
a flux CD through the ring [see the inset of Fig. 4 (a)] . Note that CD here 
only couples to the doped holes in the hopping term via the usual 
U( 1 ) degree of freedom of the conserved charge. It corresponds to the 
change of the boundary condition from a periodic one to anti- 
periodic one (CD = 7i) or a twisted boundary condition at a general flux 
CD for the doped holes (but not the spins in the superexchange term). 

Now we compute the ground-state energy difference 



A£ 2-holes 



(<D)=£ 2 (<5)-£ 2 (0). 



(3) 



If two holes are paired up, then when N x is larger than the pair size, 
one expects that A£g holes (CD) quickly vanishes at CD = n, i.e., the flux- 
quantization condition, because a pair of holes will contribute to a 2n 
phase change by winding around the closed loop once. By contrast, 
for a general twisted boundary condition of <I> < n, the contribution 
from a "coherent" Cooper pair of mass M* is expected to be 
A£2- hoks (<D)oc(Afc x ) 2 /(2M*) with Akx = Q /N]p 

Figures 4 (a) and 4 (b) illustrate the behaviors of A£g holes for the 
two-leg ladder case at CD = n and CD = nil, respectively. In Fig. 4 (a), 
with <t> satisfying the flux quantization condition <t> = n, A£g holes 
oscillates strongly and falls off exponentially as e~ N '^~ with X ~ 3. A 
roughly measures the pairing size of two holes as indicated by Ey 



shown in the same figure. On the other hand, in Fig. 4 (b) at <I> = nil, 
the coherent behavior of the hole pair as a whole is clearly shown, 
which follows an expected power-law decay ocN^ -2 at large N x [cf. 
the inset of Fig. 4 (b)]. Similar behavior has been also checked for 
other CD's deviating from the quantization n. 

The above results confirm the hole pairing at a length scale com- 
parable to the spin correlation length. Furthermore, the strong oscil- 
lation of A£ 2 - holes (7i) at small sample sizes in Fig. 4 (a) indicates a 
peculiar relative motion of a single hole within the bound pair. Note 
that previously a similar oscillation effect has been found 15 in 
A£j: hole (<D) =£i(0) -£i(0) for a single hole doped case, which is 
replotted in Fig. 4 (c) for comparison. 

As discussed in Ref. (15), such an oscillation in the single hole case 
[Fig. 4 (c)] is a direct manifestation of the so-called phase string effect 
hidden in the t-J model, which represents the non-perturbative 
quantum frustration introduced by hole hopping. Microscopically, 
the propagation of a single hole described by the t-J model can 
be precisely expressed by a superposition of quantum ampli- 
tudes of all the paths, each carrying a unique sign sequence, 
( + 1) x ( — 1) x ( — 1) x ■ ■ ■ as illustrated in Fig. 5 (a) known as 
the phase string 1113 . The sign ± in such a sequence keeps track of 
how the hole hops on the spin background by differentiating the 
microscopic processes of f- or j-spin exchanging with the hole at 
each step of hopping. The destructive interference of phase strings 
[Fig. 5 (a)], picked up by the hole from different paths, suppresses the 
forward scatterring and results in the localization of the hole based on 
the previous DMRG study 15 [cf. Fig. 4 (c)]. 

In the presence of two holes, if they form a bound pair, the strong 
frustration on the kinetic energy caused by phase strings can be 
effectively removed [as schematically illustrated in Fig. 5 (b) for 
the two-leg ladder]. In other words, the singular phase string effect 
in the t-J model provides a new non-BCS pairing force in favor of the 
charge pairing. Indeed, the hole pair behaves like a coherence entity 
in Fig. 4 (b) with the phase strings well cancelled out at large distance, 
which is obtained under a general twisted boundary condition with (D 

7t mod (n). As pointed out above, the residual phase string oscil- 
lation only shows up at the smaller N x in Fig. 4 (a), where the finite- 
size behavior of A£g holes exhibits an uncompensated phase string 
effect within the hole pair. 

As a matter of fact, such a phase string effect can be completely 
"turned off" if one replaces the hopping term H t in (1) by 



+ h.c. 



(4) 



where an extra spin-dependent sign a = ± is added, resulting in the 
so-called a 't-J model 15 . 

Then, for the a' t-J model, the corresponding AEQ ho[es exhibits no 
oscillation at all, which falls off in a power-law as 1 /N x with a = 2 
shown in Fig. 6 (a) for the two-leg ladder. Such a power-law behavior 
simply implies that the holes propagate coherently as individual 
Bloch quasiparticles 15 , without pairing. The binding energy is also 
given in Fig. 6 (a), which also decays in a power law fashion. In Fig. 6 
(b), £;, as a function of r//is present at different ladder lengths, which 
is extrapolated to a vanishingly small value in the thermodynamic 
limit. Therefore, even with the same spin gap/short-range spin cor- 
relation [cf. the inset of Fig. 6 (a)] as in the t-J case, by merely 
changing the hopping term to switch off the phase strings in (4), a 
strong hole binding is significantly reduced to negligibly small. It 
means that the superexchange interaction is by no means the sole 
pairing glue in a doped Mott insulator. 

Discussions 

In the literature it has been generally believed that the origin of 
Cooper pairing in a doped Mott insulator is solely due to the super- 
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Figure 6 | The hole binding in a two-leg ladder of the <r - (-/model, in which the phase string effect is switched off (see the text), (a) The energy difference 
= 71, is calculated at f = 3/for thegeometry showninthe inset ofFig. 4 (a). It decays in a power law fashion ~N X 2 , showing that a phase 
string free hopping term leads to a coherent quasiparticle behavior without binding. vanishes in the same power-law fashion. Note that spin 
correlations are still short-ranged as shown in the inset, (b) Binding energy for two holes as a function of f//at different sample sizes with the open 
boundary condition. The extrapolations to N = 00 show that the strong binding found for the even-leg ladders in the t- J model no longer exists here, and 
the binding strength becomes vanishingly weak. 



exchange interaction. What has been unveiled in the present DMRG 
study is, surprisingly, that the hole pairing is actually achieved by a 
combination effect. That is, the superexchange and hopping terms in 
the t-J model both play indispensable roles to the hole binding. On 
the one hand, the short-range spin correlation in a spin-gapped 
background is crucial to the spin RVB pairing. On the other hand, 
a strong frustration exerted on the kinetic energy of doped holes is 
also critical to force them to pair. 

We have also studied the a' t-J model. Its sole distinction com- 
pared with the t-J lies in that phase string signs are completely 
"switched off, while the amplitude for each path remains the same 
as the absolute weight in the latter 15 . Here, with the constructive 
interference contributed by all the paths, a coherent quasiparticle 



behavior is restored for the unpaired hole propagation. Note that 
the spin superexchange correlation still remains the same, with the 
holes having the same tendency to pair up to gain superexchange 
energy. Nonetheless, the strong pairing found in the t-J simply dis- 
appears, which unequivocally demonstrates that the non-BCS pair- 
ing force in the kinetic energy plays a critical role in the t-J model. 

So far what we have established is a novel pairing mechanism of a 
few doped holes in the f -/ladders, in which the sample size (N y ) along 
the y-axis is deliberately kept small (up to 4). The natural question is 
how much understanding that the present study provides is reliable 
on the real physics of the two-dimensional t-J model at a finite dop- 
ing, which is presumably relevant to the high-T c cuprate? Given the 
experimental facts that the Cooper pairing is more like a real space 



SCIENTIFIC REPORTS | 4:5419 | DOI: 1 0. 1 038/srep0541 9 



6 



one in a spin background of short-range AF correlations, the lessons 
learned in the present model study may be highly valuable, especially 
with the spin and charge dynamics artificially adjustable through the 
leg numbers and by turning on/off the phase string effect. In par- 
ticular, the phase string effect has been proven to be generally true for 
the t-J model on any bipartite lattice 13 . Therefore, the geometric 
limitation of the ladders in the present study is not expected to 
change the pairing mechanism fundamentally at a larger sample size 
(the leg number). 

But there is one caveat. Namely, with the increase of the leg num- 
ber, the spin gap in an undoped even-leg ladder should decrease 
monotonically, approaching to zero in the thermodynamic limit, 
where the AF long-range order is to be recovered. In the present 
study, however, we have found that a short-range spin correlation 
is crucial to the Cooper pairing. Thus the current results cannot be 
meaningfully extrapolated to the two-dimensional lattice at the same 
doping level. In other words, in the two-dimensional limit, a finite 
doping of holes will be needed in order to turn the gapless long- 
range-ordered AF state into a short-range- ordered paramagnet 
self-consistently. This is apparently beyond the scope of the present 
DMRG study. Nevertheless, our model study indicates that super- 
conductivity has to arise in a short-range "spin liquid" background, 
with high-T c (large pairing strength) emerging out of a "normal 
state" where the unpaired holes get most severely frustrated, which 
is indeed consistent with the cuprate superconductivity. 

Methods 

The numerical simulations in this work are performed by using the standard DMRG 
method 21 on both the r-/ladders in (1) and crr-/ladders in (4). Open boundary 
condition has been adopted in calculating the binding energies of the ladders, whose 
lengths have been extrapolated to the thermodynamic limit. Periodic boundary 
condition as well as twisted boundary conditions, realized by threading fluxes into the 
closed loop made of the ladder as shown in the inset of Fig. 4 (a), have been also used 
for probing the charge response. In the present DMRG simulations, we keep up to m 
— 5000 states in the DMRG block with around 20-40 sweeps to get converged results. 
The truncation error is of the order or less than 10~ 8 . 
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